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Abstract 

The equilibrium distribution function of a relativistic ideal gas has been derived 
to include the effect of angular momentum. The result agrees with the one obtained 
from kinetic theory, and consistent with relativistic thermodynamics. The role of 
angular momentum becomes transparent in this derivation, and the equilibrium 
distribution can be generalized to accommodate the effect of intrinsic angular mo- 
mentum. The results here are for flat spacetimes, however, same approach can be 
applied to static curved spacetimes. 
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The first law of thermodynamics is energy conservation. Energy plays the funda- 
mental role in thermodynamics because of its conservation, and temperature is defined 
corresponding to thermal energy. Provided that the conservation law causes the tempera- 
ture for energy, then other conserved quantities can result in thermodynamical parameters 
similar to the temperature. For example, chemical potential is introduced corresponding 
to the conservation of particle number. 

However, momentum, yet another conserved quantity, is rarely considered in ther- 
modynamics because to include it is trivial generalization with Gallilei transform in non- 
relativistic dynamics; energy can be decomposed into a simple sum of thermal energy 
and energy of bulk motion. However, such decomposition is not covariant in relativity, 
and thus momentum has non-trivial effects in relativistic thermodynamics. 

It has been proposed to define an inverse temperature corresponding to each com- 
ponent of momentum, and its three components form a four vector together with the 
inverse temperature of energy.'^''^ This four vector is called inverse temperature four 
vector, and with it we can construct manifestly covariant relativistic thermodynamics. 
The inverse temperature four vector is first introduced for a body with finite volumes 
in equilibrium,G3'l2'l3 and it is widely accepted in the relativistic thermodynamics of a 
continuum (see, e.g., Maartend^). 

Provided that momentum has a role similar to energy because of its conservation, 
angular momentum can be treated in the same way as a conserved quantity. Recently 
Nakamura^ has formulated relativistic thermodynamics with angular momentum based 
on this idea. Inverse temperature that corresponds to angular momentum is introduced 
as a covariant bi-vector (anti-symmetric rank two tensor) because angular momentum is 
a contravariant bi-vector in relativity. The purpose of the present paper is to derive this 
inverse temperature bi-vector from the view point of statistical mechanics. 

We seek for the equilibrium distribution by maximizing the entropy under the con- 
straints of conservation laws. The foundation of this approach may be highly conceptual 
and unsolved enigma in statistical mechanics. We may attribute its basis to the ergodic 
assumption of conventional statistical mechanics, or maximum entropy principle from the 
information theory as proposed by Jaynes-^^^ Whatever basis we believe, it is plausible 
that the equilibrium distribution is the one that maximizes the entropy. 

The relativistic equilibrium distribution of a uniform gas (Jiittner-Synge distribution) 
was first obtained by maximizing the entropy.'^''^''^ This approach is the same as the one 
in this paper, however, energy is the only constraint in this derivation as usually done 
in statistical mechanics. Here in the present paper we generalize it to accommodate the 
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conservation of angular momentum in addition to energy. 

Later, IsraeP has derived more general form of equilibrium distribution using kinetic 
theory. The equilibrium state is obtained as the distribution with which the collision 
term vanishes in the Boltzmann equation. This approach is different from the one by 
maximizing entropy; it is pure kinetic and entropy has no role to play there. The result 
includes distributions with non-uniform local temperature, and it was shown that the gas 
can be in equilibrium state when the local inverse temperature four vector satisfies the 
Killing equation. 

In the present paper we derive the same distribution by maximizing entropy as done 
by Jiittner or Synge.'^''^ The result includes the non-uniform equilibrium distributions 
obtained by Israel.!^ 

Temperature is derived from the Lagrange's coefficient to impose the constraint of 
energy conservation in maximizing the entropy. The same thing can be done with other 
kinds of conserved quantities; three coefficients are derived from three components of 
momentum, and these three form a four vector together with the coefficient of energy. 
This is called inverse temperature four vector mentioned above. When we apply the same 
procedure to angular momentum, the inverse temperature bi-vector is derived. We can 
also obtain the equilibrium velocity distribution, which is local Jiittner-Synge distribution 
moving along the Killing flow, in agreement with Israel.'^ 

The advantage of this approach is that the role of angular momentum conservation 
is clarified as a constraint on the entropy maximization. Also, we can apply the same 
approach to systems with gravitation, in other words, systems in curved spacetimes. 
What we investigate here is the equilibrium state in a flat spacetime, however, the same 
thing can be done for static curved spacetimes such as the Schwarzchild spacetime. If 
there exist conserved quantities in a given spacetime, we can calculate the equilibrium 
state as the maximum entropy state with the constraints of conservation laws. 

There is another advantage of our derivation; the effect of spins can be incorporated. 
The result in this paper is obtained with a model of a classical ideal gas, in which each 
particle is assumed to be a volume-less point mass. Therefore, the angular momentum 
belongs to orbital motion of particles in this model. We can extend this model to accom- 
modate intrinsic angular momentum in addition to the orbital angular momentum. 

What carries intrinsic angular momentum in this real world is spins of particles, and 
we need quantum physics to analyze it realistically. However, what we need for our 
purpose is the conservation of angular momentum only, and details of quantum aspect is 
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not necessary. Therefore, we just assume the particles carry intrinsic angular momentum 
and do not ask its details in our classical model. 

When we include intrinsic angular momentum, the sum of orbital and intrinsic an- 
gular momentum should be conserved in total. We can calculate the maximum entropy 
with the constraint of this total angular momentum conservation, and resulting equilib- 
rium state includes the thermodynamical effect of intrinsic angular momentum. Recently 
Muschik and Borzeszkowski'^' proposed an thermodynamical expression of entropy with 
the effect of spins. Their purpose is to derive a general form of entropy expression from 
thermodynamical identity, and specific form of spin effect is not given. Our result in the 
present paper is consistent with theirs, and gives its explicit form by maximizing entropy. 

The conventions in the present paper are the following. Roman indices for coordinates 
runs from 1 to 3 and Greek indices runs from to 3, where represents the temporal 
coordinate. Four vectors and three vectors are denoted by indices (e.g., x^,) and bold face 
(e.g. x) respectively, and Xq is sometimes denoted by t in non-covariant expressions. The 
metric tensor is f]f,v = diag{l, -1, -1, -1 } and the unit system is such that c (speed of light) - 
ks (Boltzman constant) - 1. 

Suppose an ideal gas is in a isolated rigid container in equilibrium. Let f{t, x, p) be 
a single particle distribution function, i.e., the probability density to find a particle near 
the point (x, p) = (x,,p,) in six dimensional position-momentum space at time t. Here 
we assume the particles in the gas are classical point masses without internal structure. 
It should be noted / is a frame-dependent expression because the time t is along the 
temporal coordinate axis of a specific frame. We need complicated form to express the 
particle distribution in a covariant way, therefore, we perform calculation assuming one 
specific reference frame in the the following. Once the equilibrium distribution obtained, 
the distribution function has the same value in other reference frames, although it is not 
a fully covariant expression.E^''^ 

What we wish to know is the equilibrium distribution that does not depend on time 
t, thus we simply denote /(x, p) = f{t, x, p) in the following. Note that this expression 
becomes frame dependent when the distribution is not spatially uniform as we will discuss 
later. The entropy is defined as 



where h is the constant in unit of angular momentum to make the argument of logarithmic 
function dimensionless. The integration is over the three volume V of the container, which 
is the cross section of its world tube sliced by the surface oit - constant. 




(1) 
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We try to find the equilibrium distribution by maximizing the above entropy under 
the constraints of conservation laws. It is known that the Minkowski spacetime has 
ten conserved quantities: four components of energy-momentum and six components 
of four dimensional angular momentum. Thus the gas has the following constraints of 
energy-momentum conservation. 



N f p^fdxdp=P^, (2) 
Jv 

and four dimensional angular momentum conservation, 

N I {x^p" - Xyp^)fdxdp = (3) 
Jv 

where P'^' , M''^ and N are the total momentum, angular momentum, and the number of 
particles respectively. To calculate ^ the temporal component of energy-momentum is 
given by p° = -yjm + pip' with m being the particle rest mass. 

The maximum entropy distribution with constraints (|2]) and dS]) can be calculated with 
the Lagrange's method as 

f = Cexp{-(}^,p^'-A^,xY), (4) 

where A^,v is an anti-symmetric tensor (bi-vector) (A^v = -Av^,) and C is a constant normal- 
ization factor to make 



I dx I dp f = 

Jv J-co 



(5) 



Note that both jS^, and A^,v are the result of global conservation laws, and thus they do 
not depend on x^,. They can be regarded as thermodynamical parameters similar to the 
temperature that characterize the whole system. 

The total entropy can be calculated by substituting @ into (d) as 

S = |6^,P'' + A^yM^'" - N In h^C (6) 

Comparing the above expression with the result of thermodynamics^' we understand that 
jS^ and A^,v are thermodynamical inverse temperatures. 

We can also confirm that the equilibrium distribution of (jU) agrees with the distribution 
obtained from the Boltzmann equation^' when we introduce the local inverse temperature 
four vector 0^ as 

d,{x) = ^, + A,,x\ (7) 
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The local velocity u^, and co-moving temperature (the temperature measured with a 
thermometer moving with the local gas velocity) T of the gas at x are then derived as 

1 6n{x) 

We can confirm i/^, is actually the local velocity of the gas by integrating the momentum 
in the reference frame moving with which gives the spatial momentum equal to zero. 

Now we have the equilibrium distribution dD, however, there was one problem in 
its derivation. We have started our calculation by assuming the equilibrium state /(x, p) 
does not depend on t = Xq, but / clearly depends on Xq when Aoi 0. This is not that / is 
truly time dependent, but that / is not spatially uniform and appears to be varying when 
measured with spatial coordinates moving relative to the non-uniform structure. 

The equilibrium state means the state when the system is isolated with the period long 
enough. When we say "time" it means the timelike coordinate in a certain coordinate 
system in relativity. When we choose an appropriate coordinate system we can obtain a 
time independent expression of (H). 

To obtain time independent distribution, we rotate the three dimensional coordinates 
so as to obtain A02 = A03 = 0, which is always possible in general. Then the factor that 
depends on time x*^ is exp[-Aoi(x°p^ - x^p'^)]. When we introduce pseudo-cylindrical 
coordinates (Rindler coordinates) for x° and x^ as 

x° = psir\hflT, x^ = pcoshflT., (9) 

then we can express the factor that depend on x° as 

exp[-Aoi(xV' - xY)] = exp(-A,ppp^) , (10) 

which does not depend on our new temporal coordinate x. 

The Rindler coordinates is known as to represent constant acceleration; a curve with 
fixed p is the orbit of the constant acceleration. Therefore, the distribution calculated 
above is the equilibrium distribution under the effect of constant acceleration. IT should 
be noted that the above expression does not cover the entire Minkowski spacetime. It is 
valid in so called "Rindler wedges", in other words, the region where the coordinate p 
becomes timelike; other regions do not have physical interest in general. 

Since the acceleration plays the same role as gravity, the above distribution has similar 
structure as the static equilibrium under constant gravity. The number density becomes 
smaller towards the acceleration direction, which is opposite to the equivalent gravitation. 
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This effect is similar to that of earth's gravity which makes the atmosphere thinner for 
higher altitude. When we take relativistic effect into account, not only the density but also 
the local temperature becomes higher towards the direction of gravity. This is because the 
energy is equivalent to the mass and becomes subject to the gravity in relativity. 

The same kind of effect takes place in a rotating body with its axis at rest. It is known 
in relativistic thermodynamics that a relativistic wheel has higher local temperature at the 
outer part; the centrifugal force plays the role of gravity. This effect can be explained from 
the thermodynamical view point, and consistent with our result in 

We will examine the effect of intrinsic angular momentum in the following. We have 
derived the equilibrium state of an ideal gas based on the angular momentum conserva- 
tion above. There the gas is modeled by classical particles without intrinsic structure. In 
contrast, real particles have spins which can contribute to angular momentum. Quantum 
theory is required for the realistic description of spins, however, we can examine qual- 
itative aspects by just assuming the existence of intrinsic angular momentum without 
paying attention to its internal mechanism. For example, one can imagine rotating clas- 
sical particles with small enough but finite volumes to mimic three dimensional intrinsic 
angular momentum. 

When the particles have intrinsic angular momentum, the particle distribution function 
/ must depend on it as /(x, p, m^^), where mf"' is the intrinsic angular momentum of a 
particle. The constraint of angular momentum ^ is then generalized as 

N f {xY - xY + m^f dxdp = M^" . (11) 

Jv 

Maximizing the entropy ^ under the above constraint and we arrive at 

/ = C exp ^-(^Y - A^v [xY + l^'")] ■ (12) 

From the above distribution function we obtain the total entropy of the gas as 

S = [^,F' + A,,{L^^ + ) - N In h'C] , (13) 

where and Z^'^ are the orbital and intrinsic angular momenta of the particles, respec- 
tively. The spatial density of intrinsic momentum in this reference frame is obtained by 
integrating NV~^mi-'^f over p and m^"^. 

The local four velocity of the gas is given by u^, in ([8]), therefore, the flux of intrinsic 
angular momentum is covariantly expressed as 

Nu r°° 

m^'Jdpdm^'^ (14) 

%J — oo 



LIV 
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where Vq is the volume of the container measured in its comoving frame. The expression 
for the local entropy flux is then expressed as 

Sp = Hp Inh'C + ej'; + ^{d,e, - (15) 

where rip and Tp are the particle number flux and energy-momentum tensor respectively 
This expression is confirmed by integrating it over V with dyd^i - dpQy = 2Af,v, which 



results in (|13|) . 



Before closing calculations, we briefly infer the relation of the above classical toy model 
to the real quantum spins. In the above model we treated the intrinsic angular momentum 
in a classical way, therefore, it has six continuous components m'^"' in (|TT]) . Consequently, 
the distribution function becomes continuous function of m/^'^ as /(x, p, mf" ). 

In contrast, what carries angular momentum in reality is quantum spins that have 
discrete values. For example, the kinetic properties of electrons (spin-| particles) are 
represented by the Wigner function (quantum distribution function) expressed as /,s(x, p) 
with spin indices r, s = 1,2; the spin indices represent two spin components in the four 
dimensional Minkowski space. Then the local angular momentum of spins cjp" is given 
by 

n 



"iYj Vp^^r{V)frs{^' P) ^P ' (16) 

where 



2 

s,r 



C = ^^»(/"-/'K(p) (17) 

with u'^(p) and yf"' being normalized sptnor components in momentum space and Dirac 
matrices, respectively; represents Dirac adjoint of u'' (see de Groot et. al.'^ for detail). 

If we can apply our entropy approach in this paper in a straight-forward manner, we 
may seek for the equilibrium distribution by maximizing the entropy Yjr,s f frs^frs dxd-p 
under the angular momentum constraint with spin contribution (|T6l) . However, we need 
careful examination to validate this procedure. For example, we have to consider the 
relation of entropy and the Wigner function f^s, but we do not know how the restriction 
of energy shell to derive f,.s affects the entropy. To examine the quantum effects further is 
out of the scope of the present paper, and we just point out the possibility to extend our 
result to quantum spins here. 

In summary, the role of angular momentum in relativistic statistical mechanics has 
been investigated in the present paper. The equilibrium distribution function dU has been 
derived by maximizing the entropy under the conservation laws of energy-momentum 
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and angular momentum. The distribution obtained is identical to the one obtained by 
Israel]^ however, the stand points on which equilibrium distributions are derived are 
quite different. 

Israel-' obtained the equilibrium state as the state in which the collision term of the 
Boltzmann equation vanishes. This approach is kinetic and the equilibrium distribution 
is obtained as the result of time evolution relaxes to the stationary state. 

Our approach here, in contrast, is to assume that the equilibrium state is the one with 
maximum entropy, in other words, the most probable state. In this approach, we do 
not pay attention to the detailed mechanisms of collisions. What we just assume is the 
existence of random energy-momentum/angular momentum exchange process; if it exist, 
we can expect the system relaxes to a unique equilibrium state that does not depend on 
actual collision mechanisms. 

To examine the angular momentum exchange involving spins, we need to consider 
quantum effects for detailed process of collisions. However, as stated above, we do not 
need it and our classical toy model here can provide qualitative insight for the equilibrium 
with angular momentum. 

The kinetic approach and entropy approach are complimentary and both have their 
own advantages. The advantage of the entropy approach is that the meaning of conser- 
vation laws becomes transparent. We have obtained the equilibrium distribution as the 
maximum entropy state under the constraints of energy-momentum and angular momen- 
tum conservation. The effect of angular momentum can be understood in this process, 
and we can incorporate intrinsic angular momentum, which is difficult to deal with the 
Boltzmann equation. We have obtained the entropy expressions with the contribution of 
intrinsic angular momentum, ((13)) and ((T5|) , consequently. 

The derivation of the entropy expressions here based on the assumption that the gas 
is in global equilibrium as a whole. However, the differential expression of (|T5)) may 
be regarded as the zero-th order entropy in the dissipative process as usually done in 
the irreversible thermodynamics. Actually, it agrees with the one assumed in the past 
literature when Op^ = 0. Further, we can expect (|T5l) is valid not only for ideal gasses 
but also other matters in general as a thermodynamical expression. If intrinsic angular 
momentum has the same effect as orbital angular momentum, the expression for total 
entropy must be in the form of ((T3|) , and then the local entropy ((T5|) can be derived. 

Recently a thermodynamical expression of local entropy is derived including the con- 
tribution of intrinsic angular momentum by Muschik and Borzeszkowski.^' Their pur- 
pose to obtain most general entropy expression starting from thermodynamical identity 



9 



without any ansatz. Therefore, the result does not give an explicit form of the entropy con- 
tribution. Our approach here is based on the ansatz of maximum entropy, and thus able 
to go one step further to determine the explicit form. The contribution of intrinsic angular 
momentum is expressed by a linear response matrix in Muschik and Borzeszkowski^ 



and its form agrees with our expression ^5} when the matrix may include differential 
operators. 

It should be noted, however, there may be some difficulties when we apply (|T5)) to local 
equilibrium with acceleration/gravity. We have obtained time independent expression 
(|TOl) by carefully choosing a curved temporal coordinate, which is possible for a global 
equilibrium. Such choice of coordinates may be difficult when each part of gas is in local 
equilibrium with different acceleration. 

We shall mention one possibility to experimentally or observationally confirm the 
results here before concluding the present paper. We have seen the intrinsic angular 



momentum can affect the equilibrium distribution as in (|12|) . What we have modeled 
here is an classical ideal gas with mimicked intrinsic angular momentum, however, this 
kind of effect is expected to take place in general. If this effect also exists in a photon 
gas, the distribution becomes dependent of the photon spin when the gas has finite total 
angular momentum. This means the black body radiation has dependency on the wave 
polarization if the radiation comes from a rotating body. Though we need quantum field 
calculation to analyze this effect in detail, a rough estimation can be made from (|T2l) : the 
effect takes place in the form of exp{-HO /T) where h is the Planck constant and O is the 
angular velocity of rotation. This may be an extremely small collection compared to the 
factor exp{-fico /T) in the Planck distribution in ordinary environments, however, highly 
accurate measurement might be able to detect it. 
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